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Abstract. We evaluate the short-range asymptotic behavior of Green's function for 
a Hamiltonian when its potential energy part has an inverse power singularity at 
the origin. The analytically solvable case of sharply screened Coulomb potential is 
considered firstly. For this potential the additional logarithmic singular term has been 
found in the short-range asymptote of the Green function as in the case of the pure 
Coulomb potential. The case of a short-range potential of an arbitrary form with 
inverse power singularity is treated on the basis of the integral Lippmann-Schwinger 
equation. It is shown that, if the singularity is weaker than the Coulomb one, the 
Green function has only standard singularity. For the case of r~ p singularity of the 
potential with 1 < p < 2 the additional singularity in the asymptotic behavior of the 
Green function appears. In the case of p = 1 the additional logarithmic singularity has 
the same form as in the case of the pure Coulomb potential. In the case of 1 < p < 2 
the additional singularity of the Green function has the form of the polar singularity 
r -p+i xhese results are applied for extending the zero-range potential formalism on 
Hamiltonians with singular potentials. 
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1. Introduction and formulation of the problem 

In this paper we concentrate on those attributes of Green's function G(z) of a 
Hamiltonian H, which are needed for evaluating the short-range asymptotic behavior 
of the Green function in the configuration space. As an operator G(z) is defined as 
G(z) = (H — z)^ 1 for z G C. The Hamiltonian is assumed to have the form 

H = -A + V(r), (1) 

where A stands for the Laplacian over r G M 3 . While the three-dimensional 
configuration space is considered, the case of arbitrary dimension d > 1 can be treated 
by similar methods. 

The potential energy term V(r) is supposed to represent a short-range potential, 
i.e. it is a real- valued smooth function for all r = \r\ > and it vanishes asymptotically 
as V(r) oc r when r — > oo. More precisely, we assume that there exists a 

constant C > such that the inequality 

\V(r)\ < C(l+r)- 1 - 6 ,5> 1 (2) 

holds for all r except of a small neighborhood of the origin r = 0. The specific concern 
of the paper is the singular short-range behavior of the potential V(r) oc r _p , when 
r — > 0. More precisely, we assume that the potential V(r) in a neighborhood of the 
origin r = can be represented as 

V(r) = r- p W(r) } (3) 

where W(r) is a smooth bounded function with the finite limit 

lim W(r) = V . (4) 

r— s-0 

In what follows this class of potentials will be referred to as the 5J(p, 5) class. For 
Hermiticity of H it is sufficient to require p < 2 and therefore it will be assumed 
throughout the paper that this inequality is fulfilled. 

The short-range asymptote of the Green function plays the decisive role in the zero- 
range potential formalism [TJ. In our recent paper [2] by studying the Green function we 
have shown that the zero-range potential has to be modified from the standard form, if 
it is constructed for the particles interacting by the Coulomb potential V c (r) = Vor -1 . 
This modification is represented in the short-range behavior of the solution <f> of the 
Schrodinger equation 
ot 

[l/r + Vblog(r)]+/3, r -+ (5) 

47T 

by the logarithmic singularity that is additional to the standard r~ l one. The 
modification appears in fact as the result of the interplay between singularities of the 
Coulomb potential and the zero-range potential. Before treating the general case we 
will deal with the sharply screened Coulomb potential V^{r) = V c (r)9(R — r), where 

| This condition for 6 can be weakened up to 6 > 0. However we will use the stronger condition S > 1 
since it guaranties the absolute convergence of integrals we deal with in the paper below. 
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R > is a screening radius and 6- function is defined as 9(t) = 1 (0) when t > (< 0). 
We consider this particular case of the screened Coulomb potential V^, which is also 
analytically solvable [3] as V c , in order to emphasize that only the short-range behavior 
of the Coulomb potential is responsible for the effect of that interplay and therefore the 
long-range behavior of the tail of the Coulomb potential does not affect the zero-range 
potential structure. 

One of the approaches for constructing the zero-range potential in the general case 
of V G 5J(p, 5) consists in inserting a delta-functional term into the Schrodinger equation 

m 

[—A + V(r) - k 2 ] 0(r, k) + X8(r)/3 = 0, (6) 

where A is a coupling constant and /3 is actually a linear functional of <fi [lj. Then the 
solution of ([6]) can be given by the Lippmann-Schwinger integral equation 

0( r , k) = Mr, k) - A / dr' G + (r, r', k 2 )5(r')(3. (7) 

Here by G + {k 2 ) we denote lim e _> G{k 2 + ie) and this notation will be used systematically 
throughout the paper. In (GO) the function (f> is the solution to the equation 

[-A + V(r)-k 2 ]Mr,k)=0, (8) 

obeying the asymptotic boundary condition 

4>o(r, k) ~ exp (ik ■ r) + Ar~ l exp (ikr) (9) 

as r — > oo. The dot-product here and further means the scalar product of vectors in M 3 . 
The integration in ((Tj) is performed easily due to the delta-function which yields 

0(r, fe) = o (r, k) - XG + (r, 0, k 2 )(3. (10) 

As it will be shown below the limit of 0o as r — > is finite for V 6 2J(p, 8) and therefore 
the non trivial short-range asymptotic of (f> is completely determined by Green's function 
term in ( jTUl) . Hence, from ( jTUl) it is seen that the principal features of the zero-range 
potential formalism follow from the short-range behavior in r of the Green function 
G + (r, 0, k 2 ) as it takes place in the case of regular potentials p]. In consecutive sections 
we will evaluate the respective asymptote of the Green function and as the result the 
zero-range potential will be constructed. 

The paper is organized as follows. In section 2 we derive the closed form 
representation for the Green function of the Hamiltonian with sharply screened Coulomb 
potential and infer the asymptote from it. As to the best of our knowledge, this 
representation for r,r' < R has been obtained here for the first time. The section 3 is 
devoted to studying the general case of potentials from 5). In the section 4 the 
asymptotes of the Green function are used for evaluating the short-range behavior of 
the solution of ([6]) and establishing the zero-range potentials for different values of p. 
Also in this section the respective pseudo-potentials are constructed. The section 5 gives 
concluding remarks. 
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2. Green's function for Screened Coulomb potential 

The Green function is defined by the solution to the inhomogeneous equation 

[-A + V£(r) - k 2 } G+{r, r', k 2 ) = 5(r - r'). 

One of the convenient ways for constructing the Green function for radial potentials is 
the use of the partial wave decomposition [4J. The Green function is represented then 
as the series in terms of Legendre polynomials Pg 

k 2 )=if:(2i+i) Gm{ry > k2) 



-PAr- r') 



£=0 



(12) 



where r = rr 1 . The partial Green function Gr£ obviously obeys the one- dimensional 
equation 



dr s 



+ 



+ V C nir) 



k 



Gne(r,r', k 2 ) = S(r - r') 



(13) 



with natural boundary condition G^i = as r = 0. The radiation boundary condition 
as r — > oo requires the outgoing wave asymptote Gju oc exp(ifcr — mt/2). 

The Green function Gri can be constructed following the standard procedure [5] 

^( r <)^( r >) 



G m (r,r',k 2 ) 



(14) 



where r- 



W(u t ,v t ) 

max{r, r'}, r< = min{r, r'} and W(v,£, Vg) means the Wronskian of solutions 



to the equation 



+ 



+ V R (r) - k 



u(r) = 0. 



(15) 



The particular solutions ue and Vg should be defined by boundary conditions ut(0) = 0, 
and vp(r) — > exp{ikr — ini/2} as r — > oo. The exact representations for both ui and 
vg depend on whether the coordinate r is in or out the interval < r < R. These 
representations can be obtained by the matching technique as in [5]. For ug one gets 

Fi(r), kr), 

(16) 



ug(r) 



r <R, 

ug{r) = aiji(kr) + bihi(kr), r > R. 
For vg the solution takes the form 

Vg{r) = aiFg{r\, kr) + b2Gg(r), kr), r < R 

Ve (r) = hj(kr), r > R. (17) 

In these representations the Zommerfeld parameter rj is defined by the standard 
expression r] = V /(2k). By ji, hi and hj we denote the Riccati-Bessel, Riccati- 
Neumann and Riccati-Hankel functions which are related to the respective spherical 
Bessel functions as for example ji(z) = zji(z). For spherical Bessel functions and for 
the regular Fg and irregular Gg Coulomb functions we use the normalization of [6]. 
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Coefficients ai, bi, a 2 and b 2 should be determined in such a way that both functions 
ug, Ve and their first derivatives are continuous at r = R. This yields 

ai = -W R (Ft, h t )/k, h = W R (Ft,j t )/k, 

a 2 = -W R (h+, G t )/k, b 2 = W R (hj } F t )/k, (18) 

where W R stands for Wronskian that is calculated at r = R. Now the Wronskian 
W(u£, Vi) from (THj) can easily be computed and takes the form 

W(ui,v i ) = W R (F l ,ht). (19) 

Equations ( |T4|) - ( TT91) completely determine the partial Green function G R e- 

For our needs of evaluating the short-range asymptotic behavior of the Green 
function G R , if R is well separated from zero, the region should be considered where 
r < R and r' < R . In this case, by inserting into f lT^|) the quantities calculated above 
we finally represent the partial Green function by the sum of two terms 

G t (r, r', k 2 ) = jF e ( v , k r< )H+( v , k r> ) + *™^F e (T), kr)F e (r,, kr'), (20) 
k k 

where XmiM) is given by 

Xm{k) = -W R (h+, H+)/W R (h+, F t ). (21) 

Here the Coulomb outgoing wave is introduced according to the definition = 
G e + iF t . 

Now we have all components which are needed for calculating the Green function 
G^ by the formula f[T2]) . In view of ( 12"U|) the representation for G^ in the region where 
r,r' < R is given by the sum of two terms 

G+{r, r', k 2 ) = G c (r, r', k 2 + iO) + Q R (r, r', k 2 ). (22) 

The first term is the conventional Coulomb Green function which is calculated by the 
partial wave series [1] as 

Gc (r, r>, e + iO) = -L f {2 ( + ^ivM^kr,) 

Airk „ — rr 
e=o 

The second term Q R reads 

Q R (r, v, *») = -L f (2£ + 1)xa(t) ^*0^faO P<(< , . (24) 

1=0 

The last formula can be rewritten in terms of the Coulomb Green functions taken on the 
upper and lower rims of the cut along the positive real semi axis of the energy complex 
plane. This can be made by using the formula Fg, = (H^ — H^)/(2i) and the method 
of the paper [3] . The result reads 



Q R (r, r') = i J dC Z*(e, C) [Gc(r, r', (, k 2 + iO)-G c (r, r', (, k 2 -ii 



1.(25) 
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Here the parameter £ is defined as £ = rr' . The kernel Zr is given by the decomposition 
in Legendre polynomials 

oo 

Zr{L = J> + mXm(k)P e (OPi(0- (26) 
e=o 

As it may be seen from the Hostler representation [7J, the Coulomb Green function 
Gji{r,r',z) actually depends on r, r' and the angle between vectors r and r' through 
the expression r ■ r'. We have reflected this fact in notations in the integrand of (125ft 
where ( stands for rr'. The formulae (I20H26P are valid only in that part of configuration 
space where r, r' < R. Nevertheless, this representation of the Green function in that 
region completely defines the transition operator Tr(z) = Vr — VrGr(z)Vr . Indeed, 
with ( 122]) the operator Tr takes the form 

T R (z) =V r - V r G c (z)Vr - V r Qr(z)Vr. (27) 

It is interesting to see what is the consequence of (|27|) when R — > oo. Evaluating the 
right hand side of ((20 asymptotically when kR > £(£ + 1) and kR > £(£ + 1) + rj 2 [8] 
we come to the expression for XRi(k) 

Xni(k) = vnexp(W e )/(kR) + 0(1/R 2 ), (28) 

where 0£ = kR — r] log(2kR) — tt£/2 + at and at = argT(£+ 1 + irj) is the Coulomb phase 
shift. From ( )26l) we have for the L,2(— 1, 1) norm of the kernel Zr 

\\Z R \\=max\xm(k)\= V /(kR) + 0(R- 2 ). (29) 

Hence, the last term in (1271) is negligible when R — > oo and therefore 

T R {z) = V R - VrG c (z)V r + 0{R- 1 ). (30) 

This can be used for calculating the limit of Tr when R — > oo and it will be done in 
another publication. 

In the last part of this section we calculate the asymptote of G + (r,0, k 2 ) when 
r — > 0. We start from the second term in (I22j) . Since [8] 

F*(77, x) = C e ( V y +1 (1 + 77x/(£ + 1) + 0(x 2 )) , (31) 
as x — > 0, the leading order behavior of Qr reads 

Q R (r, 0, A; 2 ) = C (r/)x™(fc)F (r ? , kr)/{^r) + O(fer) (32) 
as kr — > 0. From this it is seen that Q_r(^, 0, k 2 ) has the finite limit 

\imQ R {r,0,k 2 ) = kC^XRoik)/^), (33) 

where Cq{t}) = 27cr](e 27TV — l) -1 . This analysis shows us that the singular behavior 
of G~^(r,0,k 2 ) at small r comes exclusively from the first term in ( f22|) . i.e. from the 
Coulomb Green function, which as r — > has the asymptote [2J 

G c (r,0,fc 2 + i0) = i- [l/r + V \ogr] + C(k) +0(r\nr). (34) 

47T 
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Here C(k) is given by 

°W = ^ + B t log( ~ 2iA;) + V- (1 + + 2 7o - 1] , (35) 

where 70 is the Euler-Mascheroni constant and ip(z) is the digamma function [8]. 

Collecting together the expressions obtained above we arrive at the following result 
on the short-range behavior of the Green function G\ 

G+(r,0,k 2 ) = j- [1/r + Vblogr] + C(k) + ^oMXfloW + ( r ln r).(36) 

" 47T 47T 

The latter finalizes our study of Green's function for the sharply screened Coulomb 
potential. It is apparent that in this case the zero-range potential will be identical to 
that of the case of the pure Coulomb potential [2] . 

3. Green's function behavior in the case of 2J(p, S) class potentials 

In this section we study the Green function short-range asymptote for potentials of the 
2J(p, 5) class. The Lippmann-Schwinger integral equation 

G + (r, r', k 2 ) = G+(r, r', k 2 ) - J dq G+(r, q, k 2 )V{q)G + {q, r', k 2 ) (37) 



with 



Att \r — r' 

in this case has the unique solution jH |9] which completely defines the Green function 
G + {k 2 ) . For our purpose of evaluating the short-range behavior we set r' — and 
iterate (1371) one time which results 

G + (r,0,k 2 ) = G+(r,0,k 2 )- J dqG+(r,q,k 2 )V(q)G+(q,0,k 2 ) 

+ j dqG+(r,q,k 2 )V(q) j dq'G+(q, <?', k 2 ) V(q')G + (q', 0, k 2 ). (39) 

Now we consecutively consider the short-range behavior of right hand side terms. The 
first term asymptote as r — » is obvious 

G+(rAk 2 ) = ^ rr + ^ + 0(r). (40) 

For evaluating the second term on the right hand side of fl39l) it is useful to split the 
integral into two parts in order to separate the short-range and long-range contributions 
of the integrand. Let us consider the integrals 

J,(r)= f dqG+(r,q,k 2 )V(q)G+(q,0,k 2 ) (41) 

over domains flj e M 3 defined as f2i(2) = {q '■ Q. < (>)r }. The radius r can be chosen 
as any positive bounded number well separated from zero and will be specified below. 
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For the integral h{r) we can expand the Green function factors of the integrand 
by using the Taylor decomposition up to quadratic terms as 

G+(r, q } k 2 ) = l/(4vr|r - q\) + ik/(An) + 0(\r - q\ 2 ) (42) 

and similar expression holds for G^q, 0, k 2 ) with r set to 0. For the potential V(q) in 
fii we assume that r is chosen such that the formula ([3]) can be used and the W(q) 
factor can be given by its Taylor decomposition 

W(q) = V + q-VW(0) + O(q 2 ). (43) 

Then the most singular term of h{r) as r — > will be generated by introducing in (14~TI) 
with j = 1 the leading terms of integrand constituents which are defined by ( l4*2i) and 
fT43|) . This will lead us to the integral 

ij(r) = V /(4n) 2 [ dq\r- q\~ l q^' 1 . (44) 

Since we need to find the short-range behavior of this integral when r — > we can always 
take r such that r < r^. In this case the evaluation of the integral in (jHj) is easy to 
perform with the help of the formula 

i— ^-E^w-tf)' ( 45 ) 

where as usual q > = max{g, q'} and g< = min{g, q'}. The results are naturally separated 
in two ones, i.e. if p ^ 1 it reads 

'■<"> - 4,(2 \ -lf ' +1 + 4^ r/+1 <46) 
and when p = 1 the integral I± takes the form 

/?(r) = -^log(r) + ^[1 + log(ro)]. (47) 

From P6|) it is seen that if 1 < p < 2 then has the polar singularity r~ p+1 whereas 
if p < 1 the first term in (141)]) is vanishing as r — > and is regular and has a finite 
limit. From this analysis of the integral (JHJ) it becomes clear that taking into account 
the less singular terms in the expressions for the integrand of the integral I\(r) one will 
obtain non singular contributions as r — > 0. 

Let us now consider the integral l2{ r )- For the modulus of l2{ r ) we can easily 
arrive at the inequality 

| Mr) |< * /d,4^. (48) 
(47r) 2 y n2 q\r-q\ 

Let us now suppose that r is chosen such that the inequality ([2]) can be used for q > r Q 
then with the help of (j4"5|) the right hand side of fT4H|) can be estimated as 

f dq \pl)L < C [°° dqil + q)- 1 - 6 . (49) 
Jn 2 q\r-q\ J ro 

Since the last integral converges, the integral hif) is uniformly bounded for all r such 
that r < r . 
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It remains to estimate the last third term in ( I39p . The inner integral over q' by 
its structure is quite similar to the integrals considered above if G^q', 0, k 2 ) stands 
instead of G + (q',0,k 2 ). In this case the inner integral as the function of q may have 
a singularity that is not stronger than q~ p+1 one. As it has been already shown such 
a singularity in the integrand of the outer integral over q will lead to the nonsingular 
behavior of the result as the function of r in the vicinity of the point r — 0. On using 
the iterative arguments this result can easily be extended on the case of the genuine 
integrand in the third term of ( |39i) |9j. Thus, the last term in ( 1391) should have a finite 
limit as r — > 0. 

Collecting the results obtained in this section we formulate the final statement 
about the short-range behavior of the function G + (r, 0, k 2 ) for the case of 1 < p < 2 as 

G + (r, 0, k 2 ) = ±- [1/r + A)/^ 1 ] + B x + o(l), (50) 
for the case of p = 1 as 

G + (r, 0, k 2 ) = -L [1/r + V log(r)] + £? 2 + o(l), (51) 
and for the case of p < 1 as 

G + (r,0,A; 2 ) = -^ + J B3 + o(l). (52) 
47rr 

Here the constant Aq is given by 

A 



(2-p)(l-p) 

and all finite contributions from respective integrals are denoted by B~, j = 1, 2, 3. 



4. Application to zero-range potential formalism 

The zero-range potential is introduced by implementation of a special singular boundary 
condition on the solution of the Schrodinger equation at small inter-particle distances 
PQ. This boundary condition can be enforced [2] if the delta functional term is inserted 
into equation d^J). In this case the solution is represented by 

0(r, k) = <f) (r, k) - XG + (r, 0, k 2 )f3. (53) 

Here 0o is defined according to (El EJ). Since the asymptotic behavior of Green's function 
has been studied in details in preceding section it remains to estimate the short range 
behavior of 0o- It can be done on the basis of the subsequent Lippmann-Schwinger 
integral equation 

cf>o(r, k) = exp(ifc -r)-Jdq G+(r, q, k 2 )V {q)<j> Q {q, k). (54) 

It is easy to see that for a potential V G QJ(p, 5) with p < 2, 5 > 1 any iteration of ( 1541) 
is bounded function with a finite limit at r = and so is the solution. For completeness 
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let us prove this statement for the first iteration which we represent as the sum of two 
integrals J\ and J 2 defined by 

J j (r)= I dqG+(r,q,k 2 )V(q)exp(iq-k), j = 1,2. (55) 

As in the previous section, we can obtain the following estimate for the modulus of </2( r ) 
for all r such that r < tq 

Wr)|<^ f dqq- l (l + q)- l - s . (56) 

The integral on the right hand side of ( 156]) converges. Thus, ^(f) is uniformly bounded 
in the neighborhood of the origin. The integral J\{r) can be treated by similar method 
which was used in the previous section for the integral Ii(r). Particularly, if the most 
singular terms in the integrand are kept then the leading order of the integral Ji(r) 
takes the form 

Jr(r)~irf dq\r-q\- 1 q- p . (57) 



Arc 



The integral on the right hand side of this formula coincides with the integral from ([HI) 
if p + 1 in (jUj) is replaced by p. Then from the analysis made in the previous section it 
follows that the integral Ji(r) is regular as r — > 0. As the result we can conclude that 
for the potentials of the class 5J(p, S) the short-range asymptote of the wave function 
4>o(r, ft) is regular as r — > and <fto(r, k) has a finite limit at the origin. 

The singular behavior of the full solution <fi is therefore completely determined by 
the singularities of the Green function. Accordingly, we have three following cases of 
wave function asymptotes depending on the value of p. If 2 > p > 1 the asymptote is 
of the form 

cf>(r, ft) = ^ [l/r + Ao/r^ 1 } + (3, + o(l). (58) 

47T 

For the case of p = 1 it is given by 

4>(r, ft) = ^ [l/r + V log(r)] + fa + o(l), (59) 
and for the case of p < 1 the asymptote reads 

0(r,fc) = -^ + /3 3 + o(l). (60) 
Here the constants are given by following expressions 

V 



Ar 



(2-p)(l-p) 



-\(3j and j3j = O (O, ft) — \Bjf3j for j = 1,2,3. The case of p = 1 completely 
coincides with the case of the zero-range potential for the Coulomb potential [2]. 

Alternatively, the zero range potential is defined by the pseudo-potential. Following 
the procedure described in the paper [2] the form of the pseudo-potential is calculated 
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from the asymptotic expansions (EHUD])- For all three cases the pseudo-potential may 
be given by the uniform representation 

\W s {r) = A5(r)^-^ (61) 

with variables Uj determined by formulae 

= r- l + A r~ p+ \ (62) 

l/u 2 = r- l + V \og(r), (63) 

l/w 3 = r-\ (64) 



5. Conclusion 



We have shown that the Green function of a Hamiltonian with the potential possessing 
singular behavior 0(r~ p ) at small inter-particle distances receives the additional 
singularity 0(r~ p+l ) except for the case of p = 1 when the logarithmic singularity 
appears. We have considered the class 5) of potentials with p < 2 and 5 > 1. This 
choice of parameters is not exhausting. With little effort the weaker condition 5 > 
can easily be implemented in the theory by performing more delicate estimations for 
integrals I2 and J2 which should be treated as non absolutely convergent integrals. We 
have left out such an analysis in order not to make the paper too long. 
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